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I 

Abstract. From transportation networks to complex infrastructures, and to social and eco¬ 
nomic networks, a large variety of systems can be described in terms of multiplex networks 
formed by a set of nodes interacting through different network layers. Network robustness, 
as one of the most successful application areas of complex networks, has also attracted great 
interest in both theoretical and empirical researches. However, the vast majority of existing 
researches mainly focus on the robustness of single-layer networks an interdependent networks, 
how multiplex networks respond to potential attack is still short of further exploration. Here we 
study the robustness of multiplex networks under two attack strategies: layer node-based ran¬ 
dom attack and layer node-based targeted attack. A theoretical analysis framework is proposed 
to calculate the critical threshold and the size of giant component of multiplex networks when 
a fraction of layer nodes are removed randomly or intentionally. Via numerous simulations, 
it is unveiled that the theoretical method can accurately predict the threshold and the size of 
giant component, irrespective of attack strategies. Moreover, we also compare the robustness 
of multiplex networks under multiplex node-based attack and layer node-based attack, and find 
that layer node-based attack makes multiplex networks more vulnerable, regardless of average 
degree and underlying topology. Our finding may shed new light on the protection of multiplex 
networks. 
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§1 Introduction 

Robustness of networks refers to the ability of preserving their functional integration when they 

are subject to failures or attacks m- Understanding the robustness of networks is thus useful 

E-mail address: zhaodw@sdas.org (Dawei Zhao) 


1 



2 


for evaluating the resilience of systems and constructing more efficient architectures. During 
the past decades, there have been a great number of works contributing to this topic. But 
the majority of these achievements mainly focus on the vulnerability of single-layer networks 
[3Hz], which seems inconsistent with the well-recognized fact that nodes can simultaneously 
be the elements of more than one network in most, yet not all, natural and social systems 
[SHio]. Recently, Buldyrev et al. studied the robustness of interdependent networks, where 
two networks were coupled in one-to-one interdependence way Following the failure of 

one node, a cascading crash took place in both networks (namely, interdependent networks are 
intrinsically more fragile than traditional single-layer networks), which was accurately validated 
by the theoretical analysis as well. After this interesting finding, the research of network science 
is fast extended to multilayer framework [12HI6] , where systems are usually composed of several 
network layers, including interdependent networks [T7H24] . interconnected networks [25H30| and 
multiplex networks (SUISSHll]. Thus far, the topological characteristics of multilayer networks 
and dynamical process (such as evolutionary game theory |20I22) , disease spreading |28I29I37B5] , 
random diffusion m and synchronization [39) 1 upon them have attracted great attention in both 
theoretical and empirical areas (for a recent review see [12)i. 



(b) Layer-1 layer node 



Layer-2 layer node 


Figure 1: (a) six nodes are connected via two kinds of links, blue link and black link, (b) Such 
systems can be embedded into the framework of multiplex networks with two types of links. 
Each link type defines a network layer, and the nodes of each network layer are same. The 
connectivity inter-layer (dash line) is from each node to itself. 

Different from interdependent netowrks, multiplex networks, as a typical kind of topology 
structures, can be regarded as the combination of several network layers which contain the same 
nodes yet different intra-layer connections. In this sense, many real-world systems like online 
social networks [45], technological networks [46], transportation networks [47] can be further 
studied with the viewpoint of multiplex networks. Fig.l gives an illustration of multiplex frame¬ 
work: six people are connected via two kinds of relationship, for example Facebook friends (blue 
links) and Twitter friends (black links) (panel (a)). Such systems can be well embedded into the 
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framework of multiplex networks with two types of links. Each link type in the system defines 
a network layer, and the nodes of each network layer are same (see panel (b), the connectivity 
inter-layer is from each node to itself). To distinguish the node of multiplex networks (nodes in 
panel (a)) and its agent in each network layer, here we define them the terminology respectively: 
multiplex node and layer node. 

Looking back to the early topic, the research of robustness of multiplex networks thus be¬ 
comes a very interesting and crucial challenge. In [44], Min et al. explored the robustness of 
multiplex networks when multiplex nodes were removed randomly or intentionally (here the 
removal of a multiplex node means all its agents in each network layers are removed). They 
showed that correlated coupling would affect the structural robustness of multiplex networks 
in diverse fashion. In some realistic cases, however, the failure unites or attack targets may be 
just the layer nodes. For example, the users of social networks are banned to use one or some 
but not all of the social network sites. Similarly, for multiplex transport networks where nodes 
are cities and network layers are airplane network, highway network and railway network, the 
failures may take place in one or some but not all layers. Therefore, an interesting question 
naturally poses itself, which we aim to address in this letter. Namely, how does the removal of 
layer node affect the robustness of multiplex networks? 

Aiming to answer this issue, we consider the robustness of multiplex networks under layer 
node-based attack, which can be further divided into random and targeted scenarios. With the 
framework of generating function method [IS], we propose theoretical method to calculate the 
critical threshold of network crash and the size of giant component when a fraction of layer 
nodes are removed. Furthermore, we also compare the robustness of multiplex networks under 
multiplex nodes-based attack and layer node-based attack. 

§2 Model and analysis 

As mentioned in previous literatures ummu, the robustness of networks is usually evaluated 
by one critical threshold value and the size of giant component after the removal of nodes. 
If the fraction of removed nodes exceeds this critical threshold, the giant component becomes 
null. Here it is worth mentioning that the component of multiplex network is defined as a 
set of connected multiplex nodes. A pair of multiplex nodes is regarded to have connection if 
there exists at least one type of link between them. Therefore, attacking some layer nodes may 
not destroy their connection with other nodes (see Fig.2). In the following, we will focus on 
theoretical method of calculating the critical threshold value and the size of giant component of 
the multiplex networks under layer node-based attack. 

For a multiplex network composing of N multiplex nodes and m network layers, the gener¬ 
ating function for the joint degree distribution p{kj), where kj = {kj^, kj .^,..., kj^) denotes the 
degrees of a multiplex node j in each layer, can be written in the form of a finite polynomial 

m 

Go(-^)=^p(^)n4^s 


( 1 ) 
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Q Functional node 
^ Initial attacked and failure node 
Q Failure node 


Figure 2: Layer node-based attack: (a) layer node 1 of Layer-1 and layer node 4 of Layer-2 
are initially attacked, (b) soon layer nodes 1,2 of Layer-1 and 3,4 and 5 of Layer-2 become 
failure nodes since they do not belong to the giant component of corresponding layers, (c) But 
multiplex nodes 1-6 still belong to the giant component of the multiplex networks since they 
connect to the giant component through at least one type of links. 


where = {xi,X2, ...jXm) represents the auxiliary variable coupled to kj. Then the generating 
function of remaining joint degree distribution by following a randomly chosen link of network 
layer i is given by 

G«(l^) = iAGo(^), (2) 

Zi OXi 

where Zi is the average degree of layer i. 

If Mi (* = 1, 2,..., m) is defined as the probability that a multiplex node reached by following 
a random chosen link of network layer i does not belong to the giant component, it can be 
derived by the coupled self-consistency equation 

M, = GfH^), (3) 

where it = {ui,U 2 , Furthermore, the size of the giant component can be calculated 

according to 


R = 1 — 


( 4 ) 
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Along this framework, we can now turn to the layer node-based attack on multiplex networks. 
If (l)i(kj.) is used to denote the probability that a layer node with degree is removed from 
network layer i, then the generating function of the joint degree distribution after the removal 
of layer nodes can be expressed as 

m 

HoC^) = + (1 - (t)i{kj^))x\''). (5) 

t, *=i 

Correspondingly, the generating function of remaining joint degree distribution after the removal 
of layer nodes by following a randomly chosen link of network layer i is given by 

= ( 6 ) 

Zi OXi 

In the case of layer node removal, the probability Vi that a multiplex node reached by 
following one random chosen link of network layer i does not belong to the giant component can 
be written as 

Vi = — 

kj 

(I - + (1 - MkjJ)vs^‘)) (7) 

s^i 

Then, after the removal of nodes from layers, the size of giant component is given as follows 

R=l- (8) 


The existence of giant component under layer node-based attack requires the largest eigen¬ 
value A of the Jacobian matrix J of Eq. (7) at (1,I,...,1) to be larger than unity [44]. In this 
work, we mainly focus on multiplex networks composed of two Erdds-Renyi (ER) random [49j 
or Barabasi-Albert scale-free (SF) [50] network layers (namely, m = 2), J thus can be written 
as 


J = 


Kl /Cl 
/C 2 


(9) 


where Ki = [{kj1{l - (l>iikj.))) - {kj.{l - (l>%{kj.))))/Zi and K,i = {kj^kj^{l - (j)i[kj^)){l - 
4>2{,kj^)))/zi. The largest eigenvalue A is given by 


^ = 2[^i + K2 + ^/{i^ 


K2)2+4/Ci/C2]. 


( 10 ) 


§3 Results 

3.1 Layer node-based random attack 

For layer node-based random attack, which is characterized by random removal of layer nodes 
from network layers, there exists the removal probability (j)i{kj^ = {i = 1, 2; j = 1, 2, ...N). 
According to the above analysis, the critical threshold and the size of giant component of 
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multiplex networks under layer node-based random removal can be respectively expressed as 

(</)“, = 1 } ( 11 ) 

and 

= 1 - ( 12 ) 

where 4>i{kj^) = 

It is worth mentioning that above {4>i^,<p2^)c there is no giant component, whereas below 
((/)f^, 4 ’ 2 ^)c a giant connected cluster exists. 




Figure 3: (Color online) The size of giant component in dependence on removal probability 
and 4>2^ for layer node-based random attack. The black line indicates the theoretical 
critical threshold calculated according to Eq.(ll). The networks used are multiplex ER network 
with average degree (a) zi = Z 2 = I, (b) zi = 2, Z 2 = 3 and size N = 5000. 



Eigure 4: (Color online) Theoretical (line) and numerical (point) results of the size of giant 
component as a function of ^2 ^ when takes fixed values. The networks used are multiplex 
ER networks with average degree (a) = Z 2 = 1, (b) zi = 2, ^2 = 3 and size N = 5000. 

We start by inspecting how layer node-based random attack affects the robustness of mul¬ 
tiplex networks. Fig.3 shows the size R^^ of giant component in dependence on the removal 
probability and p 2 ^ for network layer 1 and network layer 2, respectively. Moreover, the 
black line indicates the theoretical critical threshold calculated according to Eq.(ll). It is clear 
that when the removal probability {pi^,p 2 ^) is above this black line, the size of giant com¬ 
ponent becomes negligible; whereas there exists one giant component if (</>i''^, ^2 ^) is located 
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below this black line. This implies that the theoretical critical threshold can accurately predict 
the impact of layer node-based attack on robustness of multiplex networks. To further validate 
this fact, we also compare the theoretical prediction derived from Eq.(12) and simulation results 
for the size of giant component in Fig.4. It can be observed that there is indeed good agreement 
between simulation and theoretical prediction. 


3.2 Layer node-based targeted attack 

Targeted attack, as a well-known attack strategy, usually aims to remove influential nodes, which 
can be identified by centrality measures, such as the degree centrality, eigenvector centrality, 
fc-shell centrality and betweenness centrality m- In this work, we mainly pay attention to the 
viewpoint of degree centrality. For layer node-based targeted attack, the removal probability of 
a layer node with degree kj^ is determined by its degree, and can be expressed as follows 

{ I, if kj ^ ^ kci 

/,, if kj. = kci . (13) 

0, if <C kci 

where k^ is the cutoff degree for attack on network layer f, and fi denotes the removal probability 
of node with degree kd- Consequently, the total fraction of removal nodes in network layer i is 
given by 

(14) 

^3 i 

where Piikj^) indicates the fraction of layer nodes with degree kj^ in layer i. 

Similar to Fqs. (11) and (12), we can get the critical threshold 

(<(>r,<(>2^^), = {(<(.r,</)niA = l}, (15) 

and the size of giant component 

= (16) 

where (j)i{kj^ is defined as Fq. (13), for layer node-based targeted attack on multiplex networks 
consisting of two network layers. 

In Fig.5, the color code represents the size of the giant component as a function of the 
removal probability and under layer node-based targeted attack, and the black line 
indicates the theoretical critical threshold calculated according to Fq.(15). Similar to Fig.3, 
the theoretical prediction fully agrees with the simulation results. Moreover, Fig.6 provides 
the further comparison between the theoretical prediction and simulation for the size of giant 
components, which also validates the accuracy of theoretical method. Combining with all the 
above phenomena, it is clear that the proposed theoretical framework can allow us to accurately 
calculate the critical threshold and the size of giant component under the layer node-based 
attack. 



Figure 5: (Color online) The size of giant component in dependence on removal probability 
and 02^ for layer node-based targeted attack. The black line indicates the theoretical 
critical threshold calculated according to Eq.(15). The networks used are multiplex ER networks 
with average degree (a) zi = Z 2 = 2, (b) zi = 2, Z 2 = 4 and size N = 5000. 



Eigure 6: (Color online) Theoretical (line) and numerical (point) results of the size of giant 
component as a function of when (pj'’^ takes fixed values. The networks used are multiplex 
ER networks with average degree (a) zi = Z 2 = 2, (b) zi = 2, Z 2 = 4: and size N = 5000. 

3.3 Comparison of robustness of multiplex networks 

Based on the above analysis, multiplex node-based attack proposed in [33], can be regarded 
as a special case of layer node-based attack when all the removed nodes or replicas are the 
same in each network layer. Erom the economic viewpoint, the cost of removing p fraction of 
multiplex nodes seems approximately equal to that of removing p fraction of layer nodes in each 
network layer. However, the damage of both scenarios on the multiplex networks may be greatly 
different. In this sense, it becomes very instructive to compare the robustness of multiplex 
networks under multiplex node-based attack and layer node-based attack. Eor simplicity of 
comparison, we assume that layer node-based attack means to remove the same proportion of 
layer nodes in each network layer in what follows. The removal probability correspondingly 
becomes for layer node-based random attack and p^'^ = pPp^ = p^'^ for 

layer node-based targeted attack. While for multiplex node-based attack, the total fraction of 
removal multiplex nodes under random attack and targeted attacks becomes p^^ (all of the 
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multiplex nodes are removed randomly with probability and 

^MT ^ (17) 

kj 

where p{kj) indicates the fraction of mnltiplex nodes with degree kj = {kj-^^ kj^ }, and 
is defined as the removal probability of multiplex nodes with degree kj and given by 




if kj^ 


%2 ^ 


</>“'(%)=< /, if%+%=fc. 


(18) 


0 , 


if kj^ 


+ < kc 


where kc is the cutoff degree and / denotes the removal probability of node which satisfies 

^jl ^j2 ^C- 



Figure 7: (Color online) The critical threshold of multiplex networks in dependence on the 
network average degree under multiplex node-based random attack (red dash line) and layer 
node-based random attack (black solid line). The networks used are (a) multiplex ER networks 
with average degree zi = Z 2 = z and (b) multiplex SF networks with average degree zi = Z 2 = z. 
The size of all the networks is iV = 5000. 



z z 


Figure 8: (Color online) The critical threshold of multiplex networks in dependence on the 
network average degree under multiplex node-based targeted attack (red line) and layer node- 
based targeted attack (black line). The networks used are (a) multiplex ER networks with 
average degree Zi = Z 2 = z and (b) multiplex SF networks with average degree zi = Z 2 = z. 
The size of all the networks is iV = 5000. 
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Similar to the above treatment, we still use the critical threshold as a uniform evaluation 
index for multiplex node-based attack and layer node-based attack. In fact, the larger the 
value of critical threshold, the better the robustness of multiplex networks against attack. Fig.7 
features how the the critical threshold of multiplex networks varies as a function of average 
degree under both multiplex node-based random attack (red line) and layer node-based random 
attack (black line). It is clear that the threshold of both cases rises with the increment of 
average degree, which means that multiplex networks are more robust for denser connections. 
Interestingly, another observation of utmost significance is that the threshold of multiplex node- 
based random attack is always higher than that of layer node-based random attack, irrespective 
of the average degree and underlying connection topology. This is to say, multiplex networks are 
more vulnerable under layer node-based attack, because it usually makes more multiplex nodes 
subject to attack and lose more connections with other multiplex nodes. Moreover, we can also 
obtain the similar observation for multiplex node-based targeted attack and layer node-based 
targeted attack in Fig.8, which further supports the fact that layer node-based attack brings 
larger damage to multiplex networks. Along this seminal finding, it may shed new light into the 
research of protection or immunization of empirical multiplex topology. 


§4 Summary 

To sum, we have studied the robustness of multiplex networks under layer node-based attack. 
Under this framework, the layer nodes can be removed randomly or intentionally, which corre¬ 
sponds to layer node-based random attack or layer node-based targeted attack. A theoretical 
method is proposed to evaluate the robustness of multiplex networks when a fraction of layer 
nodes are removed. Through numerous simulations, this method can accurately calculate the 
threshold and size of giant component, irrespective of the removal case. In addition, we also 
compare the robustness of multiplex networks under multiplex node-based attack and layer 
node-based attack. An interesting finding is that multiplex networks will be more robust un¬ 
der multiplex node-based attack, which is universal for different average degree and underlying 
topology. With regard to the reason, it may be related with the fact that layer node-based 
attack usually brings damage to more multiplex nodes, which will directly break the remaining 
joint component of networks. 

Since multiplex framework is ubiquitous in realistic social and technological networks, we 
hope that the present outcomes can inspire further research of the robustness of multiplex net¬ 
works, especially combining with the novel properties of multiplex networks, like the clustering 
characteristic [23], degree-degree correlation between network layers [37]. In addition, the tar¬ 
geted attack can also be incorporated into other centrality measures, such as the eigenvector 
centrality, fc-shell centrality and betweenness centrality [5Ij . Along this line, we may get new 
understanding for the protection of multiplex network. 
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